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Abstract 

The critical behaviour of the Ising model in the absence of an 
external magnetic field can be specified either through spontaneous 
symmetry breaking (thermal criticality) or through cluster percolation 
(geometric criticality). We extend this to finite external fields for the 
case of the Potts' model, showing that a geometric analysis leads to 
the same first order /second order structure as found in thermodynamic 
studies. We calculate the Kertesz line, separating percolating and non- 
percolating regimes, both analytically and numerically for the Potts 
model in presence of an external magnetic field. 

pacs: 05.50. +q,64.60.C,75.10.H,05.70.Fh,05.10Ln 

1 Introduction 

The critical behaviour in certain spin systems, such as the Ising model, can 
be specified in two equivalent, though conceptually quite different ways. In 
the absence of an external magnetic field, decreasing the temperature leads 
eventually to the onset of spontaneous symmetry breaking and hence to the 
singular behaviour of derivatives of the partition function. On the other 
hand, the average size of clusters of like-sign spins also diverges at a certain 
temperature, i.e., there is an onset of percolation. The relation between these 
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two distinct forms of singular behaviour has been studied extensively over 
the years, and it was shown that for the Ising model on the lattice U^, with 
d > 2, implemented with a suitable cluster definition using temperature de- 
pendent bond weights, the two forms lead to the same criticality: the critical 
temperatures as well as the corresponding critical exponents coincide in 
the two formulations [3,11]. 

In the presence of an external field H, the Z2 symmetry of the Ising 
model is explicitly broken and hence there is no more thermodynamic critical 
behaviour. Geometric critical behaviour persists, however; for T < Tp{H), 
there is percolation, while for T > Tp{H), the average cluster size remains 
finite. In the T — H plane, there thus exists a line Tp{H), the so-called 
Kertesz line, separating a percolating from a non-percolating "phase" [TT] . 
Given the mentioned correct cluster definition, it starts at Tp{0) = T^, i.e., 
at the thermodynamic critical point. 

We want to show here that the equivalence of thermodynamic and geo- 
metric critical behaviour can be extended to the case H ^ 0. Since in the case 
of continuous thermodynamic transitions, such as those of the Ising model, 
the introduction of an external field excludes singular behaviour (for the case 
d = 2 this is shown analytically [11 [6]), our problem makes sense only for first 
order transitions, for which the discontinuity remains over a certain range of 
H, even though for H ^ 0, the symmetry is broken. The ideal tool for such 
a study is the g-state Potts' model on a lattice 'Z'^, with g > 3 and d > 3. 
In this case, we have a thermodynamic phase diagram of the type shown in 
Fig. [1^, with a line of first order transitions starting at Tc(0) and ending at a 
second order point Tc{Hc) [10]; the transition at this endpoint is found to be 
in the universality class of the 3-d Ising model. In terms of the energy density 
€{H) of the system (the energy per lattice volume), the phase diagram has 
the form shown in Fig. [T]d; for H = 0, the coexistence range €2 < e(0) < ei 
corresponds to the critical temperature Tc(0). The average spin m(e) as or- 
der parameter vanishes for e > ei and becomes finite for smaller e. We want 
to show that in the temperature range Tc(0) < T{H) < Tc{Hc), the corre- 
sponding Kertesz line Tp{H) (see Fig. [T]:) coincides with that of the thermal 
discontinuity and that it also leads to the same first order/second order phase 
structure. Let us begin with a conceptual discussion of the situation. 

The g-state Potts' model in the absence of an external field provides 
g + 1 phases: the disordered phase at high temperature and g degenerate 
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Figure 1: Thermodynamic and geometric phase structure for a first order 
transition 

ordered low-temperature phases. Spontaneous symmetry breaking has the 
system fall into one of these as the temperature is decreased. Turning on a 
small external field H aligns the spins in its direction and thus effectively 
removes the g — 1 "orthogonal" low-temperature phases. Hence now only 
two phases remain: the ordered low-temperature state of spins aligned in the 
direction of if, and the disordered high-temperature phase. The two are for 
T < Tc{Hc) separated by a mixed-phase coexistence regime. At the endpoint 
T = Tc{Hc), there is a continuous transition from a system in one (symmetry 
broken) ordered phase to the corresponding (symmetric) disordered phase. 
The behaviour at H = He in Fig. [lb is thus just that of the Ising model, and 
hence the endpoint transition is in its universality class. 

In the geometric formulation for H = 0, with decreasing temperature or 
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energy density there is formation of finite clusters of q different orientations; 
the clusters here are defined using the temperature-dependent F-K bond 
weights. At e(0) = ei, the Zq symmetry is spontaneously broken: for one 
of the q directions, there now are percolating clusters, and the percolation 
strength P(e) becomes finite for e < ei. However, the disordered phase also 
still forms a percolating medium (for d > 3). A further decrease of the energy 
density reduces the fraction of space in disordered state, and for e(0) < 
e2, there is no more disordered percolation. Embedded in the disordered 
phase are at all times finite clusters of a spin orientation "orthogonal" to the 
one chosen by spontaneous symmetry breaking. In our treatment, we will 
therefore divide the set of clusters into three classes: disordered, ordered in 
the direction of symmetry breaking, and ordered orthogonal to the latter. 
While for H = 0, any of the q directions could be the given orientation, for 
H 0, the external field specifies the alignment direction, making the q — 1 
sets of "orthogonal" clusters essentially irrelevant. It is for this reason that 
at the endpoint of a line of first order transitions one generally encounters 
the universality class of the Ising model. Whatever the original symmetry 
of the system was, at the endpoint there remains only the aligned and the 
disordered ground states. 

The plan of the paper is as follows. In the next section, we recall the 
cluster treatment of the Potts' model and specify our method to identify 
the different cluster types. This will be followed by an analytic study valid 
for small external fields and by numerical calculations for different q up to 
asymptotic values of H. Formal details of the analytic calculation are given 
in the appendix. 

2 The model 

We consider a finite-volume g-state Potts model on the lattice {d > 2), 
at inverse temperature (3 = 1/T and subject to an external ordering field h. 
It is defined by the Boltzmann weight 

where the spins cTj take on the values of the set {1, . . . , g}, and where the first 
product is over nearest neighbour pairs (n,n). If we want to study the be- 
haviour of clusters, in the sense of F-K clusters, we turn to the corresponding 
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Edwards-Sokal formulation [S], given by the Boltzmann weight 

ussier, ri) = n[e-%..,o + (1 - H (2) 

where the edge variables rjij belong to {0, 1}. This "site-bond" model can be 
thought of as follows. Given a certain spin configuration, one puts between 
two neighbouring sites (Xj = aj an edge or bond with the probability 1 — e~^, 
and no edge with the probability e~^; for Ci ^ cTj, no bond is present. When 
the field is infinite, all ai = 1, and we are left with a classical bond percolation 
problem, while for finite field, one has a random bond percolation model in 
the random media given by the spin configuration. 

In the presence of an external field, we find it convenient for the study 
of the Kertesz's line to consider a modified version of the Edwards-Sokal 
formulation. We have three different types of spin combination: (0) two 
adjacent spins i,j are not equal, (1) two adjacent spins i,j are equal and 
parallel to h (we denote this direction as 1), or (2) two adjacent spins i,j are 
equal but not parallel to h. Correspondingly, we "color" the edge between 
i and j in three different colors riij, where the edge variables Uij belong to 
{0, 1, 2}. The resulting Boltzmann weight becomes 

ucEs{o-,n) = JJ[e-'^(5„,^.,o + (1 - e-^)5„,^,iX(a,=a,=i) 

i 

where the characteristic function = cij = 1) is unity for = aj = 1 
(parallel spins in the direction of h) and zero otherwise, while x(c"i = '^j 7^ 1) 
is unity for parallel spins not in the direction of h and zero otherwise. The 
summation over the spin variables then leads to the following Tricolor-Edge- 
Representation 

a;TER(n) = Y[e-^^"'^'°{l - e-'^Y^"^r^+^"^,-^^ x 

Here, Si{n) and 5*2 (n) denote the number of sites that belong to edges of 
color 1 and color 2, respectively, while 6*2(71.) denotes the number of connected 
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components of the set of edges of color 2, and |A| is the number of sites of 
the lattice under consideration. 

Let us mention that such kind of graphical representation has already 
been considered for various spin models in presence of an external field [3]. 

Let pa(^ ^ j) be the probability that the site i is connected to j by a 
path of edges of color 1 . As (geometric) order parameter we will consider the 
following mass-gap (inverse correlation length) 

m{f3, h) = — lim — In lim p/^{i ^ j) (5) 

where i and j belong to some line parallel to an axis of the lattice. As (ther- 
modynamic) order parameter, we shall consider the mean energy E{/3, h) = 
— ■|^/(/3, h), where f{P, h) is the free energy of the model @. 

3 Analytic results 

Let us first have a look at the diagram of ground state configurations of 
the TER representation which are the translation invariant configurations 
maximizing the Boltzmann weight (jlj). 

For the color a = 0,1,2, let ba be the value of the Boltzmann weight 
of the ground state configuration of color a per unit site. One finds bo = 
e-/3d(g - 1 + e'^), 6i = (1 - e-f^Ye^, ^2 = (1 - e"^)"^. 

Notice that bo = bi on the line 

/?o(/i)=ln[l + (l + (g-l)e-")i/'^] (6) 

and that bo = bi = 62 at the point /3o(0). 

The diagram of ground state configurations, inferred from the values of 
the weights bo, 61, 62 is shown in Fig. [2] (in the {h, (3) plane). 

When q is large enough and h not too large, the TER representation 
dlj) can be analyzed rigorously by a perturbative approach. Namely, by 
using the standard machinery of Pirogov-Sinai theory, we will show that the 
model undergoes a thermodynamic first order phase transition in the sense 
that the mean energy (as well as the magnetization) is discontinuous at some 
f^cih) ~ Poih). We also find for these values of the parameters, that the phase 

Note that all partitions functions, and hence the free energies, of models ([T] )-(|4] ) 
coincide 
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Figure 2: Diagram of ground state configurations: All the ground state con- 
figurations coexist at (0,/3o(0)). Below Po{h), the 0-state dominates. Above 
PoW, the 1-state dominates; it coexists with the 0-state on the line PoW, 
and with the 2-state on the line /i = 0, /3 > /3o(0). 

diagram of this model reproduces the diagram of ground state configurations 
(Fig. [2]), see Appendix for more details. 

In addition, the model exhibits a geometric (first order) transition, in the 
sense that, on the same critical line, the mass gap is discontinuous. 

Theorem 1. Assume d >2, q and h such that 

c,{l + {q-l)e~'^r'/''<l (7) 

holds, where Cd is a given number (depending only on the dimension), then 
there exists a unique Pc{h) = Po{h) + 0(1 + {q — l)e~^)~^^'^'^) such that 

1. AE{(3,{h), h) = E{l3-{h), h) - E{f3+{h), h) > 

2. m{P, h)>0 for p< (3^{h) and m{p, h) = for p > p^{h). 

The proof is given in the appendix. 
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Let us recall that it has already been shown that the Potts model ([T]) 
undergoes, for q large and h small, a first order phase transition on a critical 
line [1], where both the mean energy and the magnetization are discontinu- 
ous. Since, as already mentioned, the free energies of models ([1]) and (jl]) are 
the same, this critical line coincides with the one mentioned in the theorem. 

In the absence of an external magnetic field, the statements of the theorem 
have been shown previously [131 [HI [12] • 

Condition ([Tj) restricts the range of values of parameters to which our 
rigorous analysis applies. Moreover, we do not expect thermodynamic first 
order transitions when h is sufficiently enhanced. In the next section, we 
turn to numerical study on a wider range of values. 

4 Numerical simulations 

We have implemented a generalization of the Swendsen-Wang algorithm for 
our colored Edwards-Sokal model Q. 

First, given a spin configuration, we put between any two neighbouring 
spins of the same color, an edge colored with probability (w.p.) e~^, and 
w.p. 1 — e~^, an edge colored 1 if these spins are of color 1, and colored 2 
otherwise. When two neighbouring spins disagree, the corresponding edge is 
colored 0. 

Then, starting from an edge configuration, a spin configuration is con- 
structed as follows. Isolated sites (endpoints of 0-bonds only) are colored 
1 w.p. e^/ (g — 1 + e^) and colored c G {2, q} w.p. l/{q — 1 + e^). Non- 
isolated sites are colored 1 (w.p. 1) if they are endpoints of 1-bonds and 
colored c G {2, q} w.p. l/{q — 1). 

This algorithm allows us two compute both quantities associated to spins 
configurations and those associated to edges configurations. 

The numerical results for d = 2 are presented in Fig. [31 For g < 4, we 
found a whole geometric transition line I3c{h) for which m(/9, h) > {) when 
(3 < l3c{h), and m{j3,h) = when f3 > Pc{h)- The mass gap is continuous 
at (3c{h). For (3 < Pdh), the mean cluster sizes remain finite, while for 
P > Pc{.h) the size of 1-edge clusters diverges. The energy density as well as 
the magnetization do not show any singular behavior. 

For g > 5, some critical he appears for which the geometric transition 
line Pc{h) becomes first order when h < he, i.e. m(/3, /i) > for /? < Pdh) 
and m{P,h) = for /3 > (3c{h). In addition, on this part of the line, we 
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find that the mean energy is discontinuous 0. When h > only a geometric 
transition occurs and the scenario is the same as for q < 4. Thus our numerics 
show that the geometric and thermodynamic transitions coincide up to he, 
similarly to what we got analytically but only at (very) small field (and large 
q), see Fig. [31 

Let us mention that the numerics are in accordance with the theory for 
vanishing and infinite fields: /3c(0) = ln(l + ^g) and /3c{oo) = In 2. 
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Figure 3: Pdh) for several values of q, with "first order" behavior in red, 
"second order" in blue. The first order behavior is both thermodynamic and 
geometric. The second order behavior is only geometric. 

The system size in these calculations was L = 50, d = 2. The "first 
order" part of the transition lines has been determined via Binder cumulants 
[2]. The Hoshen-Kopelman algorithm [8] was used to study cluster statistics. 
For each value of q, more than 2 x 10^ iterations were performed. Data have 
been binned in order to control errors in measurements. 

^The Swendsen-Wang algorithm allows to compute both associated order parameters 
(mass-gap and mean energy). 
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5 Concluding remarks 



For the Potts model in the presence of an external magnetic field, we have 
shown that when the Kertesz line is first order, it coincides with the usual 
thermodynamic critical line. This property holds up to some critical point 
{hc,Pc{hc)), beyond which the thermodynamic transition disappears. Such 
behavior may well appear also for a broader class of models exhibiting first 
order transition in the presence of an external field. We believe that the 
behavior at the above critical point also belongs to the universality class 
of the Ising model, as it is the case in the 3-state Potts model in three 
dimensions |lUj . 
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7 Appendix 

We first introduce the partition function of the TER representation with 
boundary conditions a G {0, 1, 2} in a box 

where the sum is over all configurations n = {njj}jjnA05 dA is the boundary 
of A (set of sites of A with a n.n. in \ A), the notation i ~ j means that 
i and j are n.n., and 

uj^(n) = (1 - e-^Y^"^r^+^"^r^y^e-^^"^r°/^e''''^'^"^"\q - 1 + e'^)^.-^"^.^ (9) 

^The reader should not be confused by the fact that ([8]), called diluted partition function 
in PS-theory, differs from the usual one by an unimportant boundary term which makes 
the expansions (fTUl) and pT|) easier to write. 
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where ^ "1") means that the site i belongs to some edge of color 1. Next, 
consider a configurationn on the envelope of A: -E'(A) = fl A 7^ 0}. 

A site i e A is called correct if for all j ^ i, Uij takes the same value, and 
called incorrect otherwise. Denote /(n) the set of incorrect sites of the con- 
figuration n. A couple T = {Supp F, n(T)} where the support of F (Supp F) 
is a maximal connected subset of I{n), and Ti{r) the restriction of n to the 
envelope of A is called contour of the configuration n (here, a set of sites is 
called connected if the graph that joins all the sites of this set at distance 
d{i,j) = maxfc=i_...^rf \ik—jk\ < 1 is connected). A couple F = {SuppF,n(F)} 
where SuppF is a connected set of sites is called contour if there exists a 
configuration n such that F is a contour of n. For a contour F, let n-p denote 
the configuration having F as unique contour, Ext F denotes the unique infi- 
nite component of 'Z'^ \ (Supp F), Int F = Z'^ \ (Ext F U Supp F), and Int^ F 
denote the set of sites of Int F corresponding to the color m G {0, 1,2} for 
the configuration np. Two contours Fi and F2 are said to be compatible if 
their union is not connected and are called external contours if furthermore 
Int Fi C ExtA F2 and Int F2 C ExtA Fi. For a family 6 = {Fi, . . . , F„}ext of 
external contours, let ExtA^ denote the intersection A fl^^^ ExtAFfc. With 
these definitions and notations, one gets the following expansion of the par- 
tition functions over families of external contours, 

n 

Za{^)= E b[^''''''^l[p{r,) n zuint^r,), (10) 

e={ri,..,r„}oxt k=i m = 0,1,2 

where p(F) = Hiesuppr ^*("r)g^^"^^"^'''^ From ([lO]), we get 

n 

^a(A)=e' E n^'^(r,), (11) 

{Ti , . . . iPtt, }comp ^ 1 



where the sum is now over families of compatible contours and the activities 
Za{T) of contours are given by ^^(F) = p(F)6a Hm^a Hinl'^P) ' 
It is easy to prove the following Peierls' estimate, 

p(F)( max 6J-|S"PPr| < e-"IS"PPr|^ (12) 

a=0,l,2 

where e~'^ = (1 + (g — l)e~^)"^/^'^. Indeed, first notice that an incor- 
rect site i is either of color 1 or of color 2. In the first case one has 
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E,~i(^n,„o + 5n,„i) = 2d, SO that uj,{nr)/h = (e^-1) implying 
ujiinr)/ max 6a < (1 + (g - l)e-^)"(^^-'^"«^'°)/^'^. 

a=0,l,2 

Thus since 1 < X]j~'t^nij,o ^ 2c/ — 1, each incorrect site of color 1 gives at 
most a contribution to the L.H.S. of (1121) . In the second case, one has 
Ei~i(^n,„o + ^n,„2) = 2d, SO that Wi{nr)/b2 = (e^^- l)"(J^.-^"».'0)/2, implying 

LUiinr)/ max ba < (q - 1 + e'')"^^^-'^-^^'°^^^'^. 

a=0,l,2 

We then use again that 1 < X]j~i^nij,o < 2d — 1 and that C2(nr) < 
^jgg^pppX(l < 5n,3,2)/2^^~^'^"'J'^ (see [12]) to obtain that each incorrect 
site of color 2 gives at most a contribution (e'^ + q — i)-i/2+i/2d ^ g-r -j-^ ^j-^g 
L.H.S. of ([I2D. 

When the assumptions of the theorem are satisfied, the Peierls' estimate 
( IT2|) provides a good control of the system by using Pirogov-Sinai theory 
[T5] . We introduce the truncated activity 



<(r) 



^a(r)ifza(r)<(coe-oi'"pp^i 

(coe~^)l ^''PP'"' otherwise. 



where Cq is a numerical constant, and we call a contour stable if Za{T) = z'^{r). 
Let Z'^{A) be the partition function obtained from (ITT]) by leaving out unsta- 
ble contours, i.e., by taking the activities ^^(r) in fflTl) . and let us introduce 
the metastable free energies f^^\P,h) = - limAtz;''(l/|A|) InZ^(A). The 
leading term of these metastable free energies equals — Infe^. The corrections 
can be expressed by free energies of contour models which can be controlled 
by convergent cluster expansions. As a standard result of Pirogov-Sinai the- 
ory, one gets that the phase diagram of the system is a small perturbation 
of the diagram of ground state configurations. Namely, there exits a unique 
point fSciO) given by the solution of f^^P, h) = /r*(/?, h) = f^'^^p, h) for 
which all contours are stable and such that ^a(A) = Z'^[A) for a = 0,1,2. 
There exists a line Pdh) given by the solution of f^^^i^jS^h) = f'^''^{l3,h) 
when h > and such that, Za{A) = Z'^{A) for a = 0, 1. For /? < Pc{h) one 
has Zo(A) = Z'q{A), and for p > p^{h) one has Zi(A) = Z[{A). For = 
and P > /3c(0), one has in addition Z2{A) = Z^i^A). 

For the color a = 0, 1, 2, denote by h) the expectation value under 

the a-boundary condition. As a consequence of the above expansions and 
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analysis, we obtain by standard Peierls' estimates that for h> 



{K,,i)\P,h) > 1 



0{e~n for P>P,{h) 
O(e-0 for f3<f3c{h) 



(13) 
(14) 



while in addition we also get for h 



0: 



(5n.,,2)'(/3,0)>l 



O(e-0 for P>Pc{0) 



(15) 



By definition of the mean energy, one has that AE = E{/3~ , h) — E{f3~^, h) 
is proportional to the difference ('5nij,o)°(/5~, h) — (5nij,o)^(/5, h), and the first 
statement of the theorem follows immediately from these properties. 

To prove the second statement, we remark that if one imposes that the 
site i is connected to j by a path made up of edges of color 1, then under 
the boundary condition 0, there exists necessarily an external contour that 
encloses both the sites i and j. As a consequence of the above analysis 
the probability of external contours T decays like (coe)""^' ^"pp'"' when the 
0-contours are stable, i.e. when Zq{A) = Z'q{A). One thus gets p^ii ^ j) < 
(Cte e"'^)'*^''' when /5 < Pc{h) from which the first statement of the theorem 
follows. On the other hand under the boundary condition 1, the probability 
that the site i is not connected to j can be bounded from above by a small 
number 0{e~'^) when Zi{A) = Z[{A). This follows also from a Peierls type 
arguments and implies that the probability that the site i is connected to j 
under the boundary condition 1 is greater than 1 — 0(6""^) for j3 > l3c{h). It 
gives also that the probability pA^i ^ j) for the site i to be connected with j 
under the boundary condition is also greater than 1 — 0(6""^) for P > Pc{h), 
implying the second statement. 
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